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A new numerical method is proposed for determining the low-frequency dynamics of the charge 
carrier coupled to the deformable quantum lattice. As an example, the polaron band structure 
is calculated for the one-dimensional Holstein model. The adiabatic limit on the lattice, which 
cannot be reached by other approaches, is investigated. In particular, an accurate description is 
obtained of the crossover between quantum small adiabatic polarons, pinned by the lattice, and 
large adiabatic polarons, moving along the continuum as classical particles. It is shown how the 
adiabatic contributions to the polaron dispersion, involving spatial correlations over multiple lattice 
sites, can be treated easily in terms of coherent states. 

The motion of charge carriers, coupled to the phonon field on a lattice, has been an object of intensive investigations 
for many years. In this context, different types of quasi-particles are studied, such as polarons, bipolarons, topological 
solitons, to name a few. In particular, polarons were first introduced as states of broken translational symmetry, 
describing a localized electron trapped by the self-induced, static lattice deformation^. By taking the lattice kinetic 
energy into account, the translational symmetry is restored: the deformation moves along the lattice followed almost 
instantaneously (adiabatically) by the electron, provided that it is fast enough. For large adiabatic polarons, which 
behave as free particles, the effective mass can be obtained in the continuous approximation, treating the lattice 
classically^. In the opposite (small polaron) limit, the polaron is strongly pinned by the discreteness of the lattice. 
For this latter case it is well-known that the quantum treatment is necessary to describe the dispersion properly^. 
Although the translationally invariant diagrammatic perturbation theory^ or similar methods provide in principle the 
unifying approach to the crossover between large and small (adiabatic) polarons, the corresponding calculations are 
intricate and have not yet been carried out. 

This long-standing problem can be simply solved by the relevant coherent states method (RCSM) developed here. 
As an illustrative example, the RCSM is applied to the ID Holstein polaron, involving an electron coupled locally to 
the lattice displacement. The Holstein Hamiltonian is given by^ 

H ^ -t'^ciiCn+l + Cn-l) + UJo'^blbn - g^cicnibi + bn) , (1) 
n n n 

with Cj (cj) and b^j (bj) the creation (annihilation) operators for the electron and phonon at the site j, respectively. 
t is the nearest-neighbor electron hopping integral, ujq is the energy of the dispersionless optical phonon branch, 
and g is the electron-phonon coupHng constant. There have been extensive studies^ of the polaron states for the 
Hamiltonian Accurate results were obtained for values of the adiabatic ratio t/wo ^ 5, which are far too small 
to address the formation of the large adiabatic polaron. On the other hand, the RCSM does not suffer from such 
restriction. With modest computational effort, it can be applied for any value of t/ujQ, including the t/ujQ ^ 1 limit. 
In this way, the crossover between small and large adiabatic polarons can be investigated. 

The RCSM is basically a variational procedure performed in two steps. In the first one, a finite set of localized 
polaron wave functions |(^|) is selected, capable of describing the polaron adiabatic motion within the unit cell 
surrounding the site j. For the site j, the localized polaron wave function is expressed as their linear combination 
= X^s^sl'/'j)' '^ith variational coefficients a^. In the second step, the translational invariance of the problem is 
restored by writing the polaron wave function in the translationally invariant form, 

with K the polaron momentum. The set of wave functions is linearly independent provided that ^ |(y9y) 

for all j 7^ j' and s ^ s' . The minimization of the variational energy with respect to the complex coefficients a* yields 
the generalized eigenvalue problem, 

s' s' 

the solutions of which are the lowest and excited polaron states for a given momentum K. 
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For a given set of parameters the RCSM trial wave function ([2]) is chosen according to the properties of the effective 
potential that characterizes the adiabatic dynamics of the polaron. This makes the RCSM essentially different from 
the other variational approaches in the literature. Each localized polaron wave function in eq. ^ is formed as a 
product of the electron and the lattice part, 



S,{x^+if)\Q) . (4) 

The wave function (0]) is parameterized by the iV-dimensional complex vector + if? = where N is the 

number of lattice sites {N oo). The lattice part in eq. ^ is given in terms of the coherent states, i.e., of the 
displaced harmonic oscillators, 

n 

Here, the vector component defines the mean lattice displacement at the site j + n, whereas defines the mean 
lattice momentum. In eq. (|4]), the electron state rjnix'^) is given by the electron ground state corresponding to the 
mean lattice deformation In this respect, the method of calculating rjnix'^) in eq. |[4]) follows the same steps as the 
adiabatic approximation^. In particular, the electron wave function is obtained as the ground state ly = solution of 
the equation 

e^-^ Vi^^ = t {r^[:l, + ri\:\) 2.9 x^ . (5) 

With rjn = rin'^ in eq. ([4]), Ylm k/iP = 1) the expectation value of the Hamiltonian ([1]) is obtained for the localized 
polaron wave function (0]) as 

{v^^\H\^^)=oj^\f\^ + ^,\x^\'+e{x^) , (6) 

where £(x^) = e^^\x'^) is the electron ground state energy of eq. ([5]). The last two terms in eq. ([6]) define the adiabatic 
potential Uad{x), i.e., the potential energy of the lattice in the adiabatic approximation. This simple relation between 
the wave function ^ and the adiabatic potential provides a well-controlled and efficient way of treating the adiabatic 
correlations. 

Although x" can represent any lattice deformation in eq. (j4]) , it is important to realize that the set of wave functions 
\(Pj) in eq. ([2]) is meant to capture the polaron dynamics around the site j, whereas the delocalization of the polaron 
over the entire lattice is accounted for by the translationally invariant form of the trial wave function ([2]). Furthermore, 
it should be stressed that the trial wave function ^ includes nonadiabatic processes as well. Namely, each localized 
polaron wave function |(^|) in eq. (0]) fulfills the adiabatic approximation only to the mean-field level. That is, in eq. (0]) 
the localized electron wave function (a;* ) is determined through eq. ([5]) by the mean lattice deformation x^ . On the 
other hand, the lattice quantum fluctuations around x'^ associated with coherent states, in eq. (0]) are not followed 
adiabatically by the electron. In particular, one can easily verify that the RCSM reproduces the well-known small 
nonadiabatic polaron^. For A = Sp/t » 1, the polaron lattice deformation localizes adiabatically to a single lattice 
site, x^ w ff/wo Sn,o, where Snfl is the Kronecker delta. By considering solely this lattice conflguration corresponding 

to the minimum of the adiabatic potential, eq. ^ gives T]n^ w Sn,o, which yields, through eq. ((2l), the polaron wave 
function as 

3 

The corresponding expectation value of the energy ([3]) is equal to the small nonadiabatic polaron energy — £p — 
£k e~^p/'^°, with Ek the free electron energy. Such dispersion of the polaron is purely nonadiabatic: during the 
hopping to the nearest-neighbor sites the electron and the lattice deformation detach from each other. 

When the adiabatic contributions to the polaron dispersion, described by a joint motion of the electron and the 
lattice deformation, are dominant the situation is fundamentally different. In particular, the adiabatic motion of the 
polaron between lattice sites is characterized by two types of stationary points of Uad{x), the minima xm and the 
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saddle points with the lowest energy xpN (i.e., two types of lattice deformations). The first correspond to polarons 
centered at the lattice site, whereas the second correspond to polarons centered half-way between the lattice sites. 
Their difference in energy defines the Peierls-Nabarro (PN) barrier Ap^r, 



ApAT = Uad{xpn) - Uad{xm) , 

which is the minimal energy barrier that must to be overcome in order to move the polaron classically and adiabatically 
from one lattice site to another. 

The adiabatic motion of the large polaron within one unit cell involves changes of the lattice deformation over many 
lattice sites. This difficulty of dealing with a large number of lattice coordinates can be avoided by introducing a 
smaller set of new coordinates that are most relevant for the polaron dispersion. One of the possible solution is to 
work with the moving set of coordinates ^ and Qq(C), where ^ is the centroid coordinate determining the position of 
the polaron along the minimal energy path of its motion, while Qa{C) are the coordinates of the normal modes a that 
are moving with the polaron^ii^^. The efficiency of this approach is derived from the fact that the effective coupling 
of i with QaiO is significant only for a few modes a. 

While in the RCSM context, the moving set of coordinates can be used quite generally for various electron-phonon 
models. An alternative choice of coordinates is considered here for the Holstein model, based on the specific properties 
of the adiabatic potential. The homogenous lattice deformation in eq. HI) couples only to the total electron density. 
Consequently, the harmonic dynamics of the homogenous q = lattice mode can be separated out, and one has 
to consider explicitly only the lattice deformations satisfying the sum rule J^n^n = g/^iy^- For these physically 
interesting deformations x — g/ujQ x, where x denotes the A^-dimensional unit vector, it follows from eq. (HJ that 
Uad{x) exhibits a simple scaling behavior, to our knowledge not noted before. 



Uad{x) = wol^P + e{x) = Ep Uad{x, A) , (7) 

with Ep = g'^/ujQ and A = Ep/t determining respectively the amplitude and the shape of the adiabatic potential. One 
immediately recognizes that A governs the unit vectors corresponding to stationary points, given by ^Uad{x) — ^ 
Xn = g/i^o hnP- These points can be found^ from the discrete nonlinear Schrodinger equation |[5|). Furthermore, 
A governs the harmonic adiabatic dynamics around the minima of Uad{x). The expansion of Uad{x) in powers of 
displacements 5 from the equilibrium deformation xm , 5 = x — xm , '^s. given by 



Uad{.x) = Uad{xm) + Wo^n„_™ dnSjn + UaE ■ 

n.m 

The second term determines the harmonic fluctuations, while the adiabatic anharmonic dynamics is described by the 
higher order contributions . The Hessian matrix depends on A, and it can be expressed by means of the 

electron basis ([5]) corresponding to xm^'^, 



with 5n,in the Kronecker delta. Up to the second order in <5, the lattice vibrations decouple in terms of normal modes 
a, which are characterized by normal coordinates 5a and eigenenergies uoa- Notice that both, the unit displacement 
5a along the normal coordinate 5a and the ratio uja/oJo, with ujq the bare phonon energy, are determined by A. 

As shown here, the use of normal coordinates 5a provides a simple and efficient way of keeping track of the 
lattice configurations that are relevant for the motion of the Holstein polaron associated with one unit cell. Namely, 
by considering the displacement from the equilibrium configuration to the nearest saddle point in terms of normal 
coordinates. 



5pN = XpN - Xm ^ g/l^O 2^ C0S7q 5a , 

a 

one finds that 5pn lies almost entirely in the plane of the two lowest normal modes. That is, for all values of parameter 
A, which determine the arguments 7q., one obtains cos^ 7p + cos^ > 0.99. Here, P (pinning) and B (breather) are 
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Figure 1: The adiabatic potential for the small polaron case (A = 2), as a function of coordinates 5p and 5b, belonging to the 
pinning and breather normal modes. The absolute minimum, the two saddle points and the two local minima are denoted by 
M, PN, and Ml, respectively. 

used to denote the two lowest, odd and even, normal modes. A better insight into the properties of the adiabatic 
potential can be obtained from fig. [H which shows the adiabatic potential for A = 2 as a function of displacements 
along the normal coordinates Sp and Sb- The absolute minimum of the adiabatic potential is at the origin, while the 
saddle points PN correspond to the top of the PN barrier. It is satisfactory that the two local minima are given by 
lattice configurations which are close to the equilibrium deformation of the polaron centered at the nearest-neighbor 
sites. 

The above findings related to the properties of the adiabatic potential imply that in order to describe the adiabatic 
polaron motion properly it is sufficient to take into account the adiabatic correlations related to the pinning and 
breather modes in eq. ^ . Such a procedure breaks the translational symmetry of the adiabatic potential ^ , 



U{Sp, 6b), shown in fig.IU describes the part of the adiabatic potential Uad{x) that is relevant for the polaron motion 
within one unit cell. The adiabatic correlations related to higher normal modes, a ^ P,B, which weakly affect the 
polaron dispersion properties, are not of significant physical interest here. These nearly harmonic correlations result 
in the softening of the higher modes. If not ignored, this softening would be manifested by a set of eigenenergies Ua 
that replace ojq in the third term of eq. ([8]), uja ^ i^o- 

In the present case, having identified the two essential coordinates for the adiabatic motion of the Holstein polaron 
within one unit cell, Sp and Sb, it is straightforward to proceed with the RCSM. One considers a set of wave functions 
(U given by complex amplitudes + ijr' corresponding to the points on the grid lying in the plane of the two lowest 
normal modes. 



with n,m,n' ,m' , integers, and Ax, Ap, grid spacings. The finite set of states (0]), important for the low-frequency 
dynamics, can be chosen by considering only such x" and for which the energy ^ is below a certain value. For 
the polaron ground state, this value should be at least of the order of a few uiq, or of the order of the PN barrier if 
this barrier is largeii. The grid spacing Ax in eq. ([9]), as it can be seen from eqs. ((2| and ((4|, defines the number 
of different electron wave functions rjn involved in calculations. The adiabatic correlations are described with better 
accuracy as this number is increased by decreasing Ax. However, when Aa; is too small. Ax ^ 1, some functions 
(|4]) are nearly identical, and eq. ([3]) cannot be solved numerically. Nevertheless, satisfactory results are obtained for 
Aa;, Ap ~ 1 already. 

In order to check the vaHdity of the RCSM, a comparison with the results of the previously developed^^ exact 
translational method (ETM) is performed. For i/wg ^ 5 the ETM provides the polaron band structure below the 
phonon threshold with great accuracy^ii^. This makes it particularly suitable for the analysis of the RCSM, which 
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x'^ — xm + n Ax Sp + m Ax Sb , jf — n' Ap Sp + m' Ap Sb , 



(9) 



Figure 2: Panel a): the polaron band structure below the phonon threshold for t = 5 resolved by the ETM (solid curves) and 
the RCSM (dashed curves). Panel b): ETM vs. RCSM polaron effective mass Jriej/mpjj for the two lowest i = 0, 1 bands 
(t = 5). Panel c): the polaron band structure for t = 250 in the crossover regime by means of 22 states with different momenta 
for each band, K = 0,-K,n x 0.15, n < 20. tJo is used as the unit of energy. 

exhibits the excited polaron bands as well. Figure [2k shows the lowest i = and the two i = 1,2 excited bands 
related to the pinning and the breather mode in terms of the band boundaries corresponding to K = and K = tt 
polaron states. All the curves are given with respect to the ETM ground state energy. The solid curves are the ETM, 
whereas the dashed are the RCSM results. The choice of parameters in fig. [2k corresponds to the crossover regime 
between free and pinned polaron states. As discussed in detail in refj^, for g/uiQ w 3.2 there is a strong hybridization 
of the i = 1 and i = 2 excited bands. Such behavior may be explained from eq. |(8]) and fig. [B the two lowest normal 
modes are effectively coupled by the polaron motion. The ETM bands shown in fig.[2k are reproduced very accurately 
by the RCSM. The differences become pronounced only close to the phonon threshold E'^^'' « ivq, where long-range 
correlations mediated by weakly correlated phonon excitations become important due to the existence of the continuum 
of states above wcr^. Figure [2}d shows the inverse effective mass of the polaron, Tiej/m^^jj = SS^'/^^^f l^=o for the 
same parameters as in fig. [2^. One finds that the agreement between the ETM (solid curves) and the RCSM (dashed 
curves) is excellent for the i = 0, 1 bands, even though nipoi varies over four orders of magnitude. 

The polaron crossover from pinned to free polaron states, shown in fig. [2k for t/ujQ ~ 5, corresponds to an adiabatic- 
nonadiabatic crossover in the polaron dynamics. Namely, all the excited bands below the phonon threshold, related 
to the soft adiabatic phonon modes, at the weak-coupling side of fig. [2k shift towards the phonon threshold, indicating 
that the adiabatic correlations become suppressed completely by the nonadiabatic contributionsi^. On the contrary, 
the polaron crossover shown in fig. [2]: for t/ujo = 250 is governed entirely by the adiabatic dynamics. Figure[2l: shows 
the lowest and the excited bands corresponding to the pinning and breather modes, shifted by the RCSM ground state 
energy. Each of the bands is given by the 22 states with different momenta, K = 0,TT,n x 0.15, n < 20. On the left 
side of fig. [21: the polarons are free (large), i.e., the PN is neghgible. The polaron dispersion is approximately given 
by Ek ~ for the whole energy range shown in fig. [2]:. For g/ujo « 15 the energy of the K = state involving the 
excited breather mode is very close to the value derived in the continuous adiabatic approximation (large adiabatic 
polarons), ujb/i^o = 0.65. On the right side of fig. [2b the spectrum is characterized by a large PN barrier (small 
adiabatic polarons). That is, the bandwidths are exponentially small, while the position of the excited bands are 
given by the multiples of the pinning and breather mode energies, mup + m'uiB- The crossover between pinned and 
free adiabatic polarons takes place in the central part of fig. [2fc, where strong hybridization between excited bands 
occurs due to the pinning effects. 

It is important to realize that the polarons at the left side of fig. [2k are fundamentally different from the polarons 
at the left side of fig. [2k. While the latter are nonadiabatic and characterized by strong quantum effects, the effective 
mass of the former is predicted correctly by the classical soliton theory, i.e., as for the classical lattice deformation 
that moves along the lattice. The band structure in fig. [2k, exhibiting many polaron bands, has never been obtained 
before. It reveals how the dispersion of the adiabatic polarons changes its character from the quantum to the classical 
one, as the PN barrier decreases. 

In conclusion, the newly proposed RCSM treats the low-frequency dynamics of charge carriers which interact with 
the quantum lattice. Its application to the ID Holstein polaron problem provides the polaron band structure in the 
adiabatic limit t ^ ujq, not reachable by any other available method. In particular, the RCSM gives an accurate 
description of the crossover between the small and large adiabatic polarons, for which the nature of the polaron 
dispersion changes from quantum to classical one. For t ^ luq the accuracy of the RCSM is verified by comparison 
to the practically exact ETM results. The only differences of note are found only for weak couplings close to the 
phonon threshold when the long-range nonadiabatic correlations are important. As the weak-coupling regime can be 
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described by the perturbation theory (PT) or ETM for any t/wo, one finds that the RCSM together with the PT or 
ETM provides a solution of the Holstein polaron problem for the entire parameter space. 
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